Fano resonances with discrete breathers 
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A theoretical study of linear wave scattering by time-periodic spatially localized excitations (discrete 
breathers (DB)) is presented. We obtain that the wave propagation is strongly influenced by a local 
coupling between an open and closed channels generated by the DB. A peculiar effect of total 
reflection occurs due to a Fano resonance when a localized state originating from closed channels 
resonates with the open channel. For the discrete nonlinear Schrodinger chain we provide with 
an analytical result for the frequency dependence of the transmission coefficient, including the 
possibility of resonant reflection. We extend the analysis to chains of weakly coupled anharmonic 
oscillators and discuss the relevance of the effect for electronic transport spectroscopy of mesoscopic 
systems. 
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It is a well established fact that various nonlinear spa- 
tially discrete systems can support time-periodic spa- 
tially localized excitations called discrete breather states 
(DB) [1]. These states originate from a peculiar inter- 
play between the nonlinearity and discreteness of the lat- 
tice rather than from a disorder. While the nonlinearity 
yields an amplitude-dependent tunability of frequencies 
of DBs, fib, the spatial discreteness of the system leads to 
finite upper bounds for the frequency spectrum of small 
amplitude plane waves oj q . This tunability allows one 
to escape resonances of all multiples of the breather fre- 
quency Qb with the plane wave frequencies u q , and cor- 
respondingly to stabilize the DB state. The frequency 
dependent localization length of DB's and their stability 
with respect to small amplitude perturbations have been 
widely studied [1]. DBs have been observed in experi- 
ments covering such diverse fields as interacting Joseph- 
son junctions [2], magnetic systems [3] and lattice dy- 
namics of crystals [4] . 

For propagating linear waves a DB acts as a time- 
periodic scattering potential, and the transmission coef- 
ficient T depends on both the wave vector q of the linear 
wave and the breather frequency fib- The most pecu- 
liar effect, observed in many numerical studies of wave 
scattering by DBs, is the total reflection as T = [5,6]. 
Note that the presence of a static potential cannot lead to 
such a total reflection in one-dimensional systems. Sim- 
ilar features are also discussed in other areas, such as 
electron transport through point contacts, quantum dots 
and wires [7,8]. The crucial condition allowing a total 
reflection in these systems is the presence of a few cou- 
pled channels connected with the transverse direction of 
motion. On the other hand, the wave propagation in the 
presence of a time-periodic scattering potential is char- 
acterized by open and closed channels emerging from the 
Floquct formalism [1,5,6]. The open channel guides the 



propagating waves, while the eigenfrequencies of closed 
channels do not match the spectrum of linear waves. 

In this Letter we show that the total reflection of linear 
waves in the open channel occurs when a localized state 
originating from one of the closed channels resonates with 
the open channel spectrum, a condition similar to the 
well-known Fano resonance [9]. We will use this un- 
derstanding to predict Fano resonance positions for wave 
scattering by breathers in weakly interacting anharmonic 
oscillator chains and to discuss the relevance of this effect 
for electronic transport spectroscopy. 

We start our study of the wave scattering by DB with 
the discrete nonlinear Schrodinger system (DNLS) that 
has been used frequently to study breather properties due 
to its tractable form. Wave scattering by breathers in the 
DNLS was studied numerically in [10], where resonant 
total reflection was also observed. First, we provide an 
analytical solution for the DNLS scattering problem. The 
equations of motion for the DNLS are 

= C(tf n+ i + *„_].) + |*„| 2 *„ (1) 

where the integer n labels the lattice sites, ^ n is a com- 
plex scalar variable and C describes the nearest neighbor 
interaction on the lattice. For small amplitude waves 
\Jr n (i) = ee «("9*-9» l ) the dispersion relation is 

LO q = — 2Ccosq . (2) 

Breather solutions have the form 

*„(*) = A n e~ inbt , i, nHoo -> (3) 

where the time-independent amplitude A n can be taken 
real valued, and the breather frequency Qb ^ Lu q is & func- 
tion of the maximum amplitude Aq. The spatial local- 
ization is given by an exponential law A n ~ e~ A '™ where 
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cosh A = |f2j,|/2C [11]. Thus the breather can be ap- 
proximated as a single-site excitation if \Clb\ >• C. Then 
the relation between the single-site amplitude A and 
becomes Q& = Aq. In the following we will neglect the 
breather amplitudes for n^O, i.e. we assume A n ^o ~ 0, 
since A±i w jf-Ao *C -A . We emphasize that in the limit 
of a nearly single-site localized breather solution in the 
DNLS the solution is linearly stable [11]. 

We add small perturbations to the breather solution 

*„(*) = *»(*) + <t>n{t) (4) 

and linearize equations (1) with respect to <j>n{t) : 

ij> n = C(<f> n+1 + </»„_!) + n b 8 nfi (2cj> + e- 2in *VS) (5) 

with 5 n ,m being the usual Kronecker symbol. The general 
solution to this problem is given by the sum of contribu- 
tions due to two channels, 



b n {t) = X n e luJt + Y^g-^e-Mt 



(6) 



where X n and Y n are complex numbers satisfying the 
following algebraic equations: 

-LoX n = C(X n+1 + X n _i) + Sl b 6 n , (2X + Y ) , (7) 
(2il b + u;)Y n = C(Y n+1 + r n _!) + fl b S nfi (2Y + X ) . (8) 

Away from the breather center n = 0, Eq. (5) allows 
for the existence of plane waves with the spectrum u> q . 
Keeping in mind the propagation of waves, we set ui = uj q 
for some value of q. Thus the X channel is an open one, 
while the Y channel is a closed one, i.e. its frequency 
— (2fi(, + uj q ) does not match the spectrum u> q . 

Instead of solving Eqs. (7) and (8) we will consider a 
more general set of equations 

-u q X n = C(X n+1 + X n _i) - 8 n , (V x X + V a Y ) , (9) 
{SI + oj g )Y n = C(Y n+1 + y n _!) - 6 n ,o(V y Y + V a X ) (10) 

which is reduced to (7) and (8), if O = 2il b and V x = 
V y = 2V a = — 20fc,. For a particular case V a = 0, i.e. 
when the closed V-channel is decoupled from the open 
one, the former possesses exactly one localized eigenstate 
due to a nonzero value of V y with the frequency 



= -n + ^v 2 + ac 2 . 



(ii) 



To compute the transmission coefficient T we use the 
transfer matrix method described, e.g. in Ref. [12]. 
The boundary conditions arc X^+i = re lq , Xjy = 
t , Y/v+i = D/k , Yjv = D for the right end of the chain 



and X-n-1 = 1 



X. 
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F , Y"_at = kF for the left one. Here r and p are the 
transmission and reflection amplitudes with T — |t| 2 = 
1 — \p\ 2 . F and D describe the exponentially decaying 
amplitudes in the closed V-channel, where the degree of 
localization is described by the coefficient k = e A : 
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2C 



n + Lo q + ^{si + u q y - Ac 2 



(12) 



The 4x4 transfer matrix is defined by Eqs. (9) and (10) 
at n = 0. After finding the solutions of the corresponding 
four linear equations we obtain 



T = 



4 sin q 



(2cosg-a-^r) +4sin 2 g 



a = 



-LU q + V x 



b = 



Q + LO q + Vy 



d = 



V a 



(13) 



(14) 
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This is the central result of this paper, which allows one 
to conclude that total reflection is obtained when the 
condition 



2 - bn = 



(15) 



is realized. It is equivalent to the resonance condition 

uj q = Jv\ (16) 

which has a clear physical meaning: total reflection oc- 
curs when a local mode, originating from the closed Y- 
channel, resonates with the plane wave spectrum uj q of 
the open A-channel. The only condition is that the cou- 
pling between the open and closed channels V a is nonzero. 
Remarkably the resonance position does not depend on 
the actual value of V a , i.e. there is no renormalization. 
The existence of local modes, which originate from the 
X-channel for nonzero V x and possibly resonate with the 
closed Y"-channel is evidently of no importance. Eq.(13) 
also yields zero transmission for q = 0, tt due to a van- 
ishing of the group velocity du> q /dq ~ sing for these q 
values; we will not focus on these trivial total reflections. 

This resonant total reflection is very similar to the 
Fano resonance [9], since it is directly related to a local 
state resonating and interacting with the continuum of 
extended states. The fact that the resonance location is 
independent of the coupling V a is due to the local, single- 
site, character of the coupling between the local mode 
(originating from the F-channel) and the open channel. 
If this coupling has a finite nonzero localization length, 
i.e. several neighboring sites are involved, then the total 
reflection happens at an energy, which does not neces- 
sarily coincide with that of the local state (16) [13]. A 
more physical formulation for the condition of absence 
of significant renormalizations of the resonance location 
is that the wave length of the propagating wave is large 
compared to the extension of the space region where the 
channel coupling occurs. 

Returning to the case of a DNLS breather at weak 
coupling, we insert the values for Q , V x ,V y and V a 
into (13,14) and obtain the following expression for the 
transmission coefficient: 



4 sin q 



( 20(, 

V c 



oT^r tt- 2 + 4 sin 2 q 

2C £ 1+K cos q I * 



(17) 
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The result is that any breather solution of the DNLS close 
to the anticontinuous limit (the interaction C goes to 
zero) provides us with a total reflection in a close vicinity 

C 

of q = 7r/2. Indeed, if we expand (17) in — we obtain 



fib 



4 ^ 4 .In 

-^-sin 2q 



c 



(18) 



in the lowest order, provided -C | cosg|. 
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FIG. 1. Transmission coefficient T for DNLS breather with 
C = 0.01 and fib = 1 versus q. Dashed line - numerically exact 
result, solid line - analytical result (17). Inset - zoom around 
total reflection, where dotted line marks q = it/2. 

Fig.l compares a numerically obtained q dependence 
of the transmission coefficient with (17) for C = 0.01 [14]. 
We obtain a very good agreement, except for a small shift 
of the true total reflection position with respect to q = 
7r/2. It is due to the small but nonzero finite extension 
of the scattering potential, which leads to a spread of the 
coupling between the local mode and the open channel. 

Next we consider the system of weakly interacting an- 
harmonic oscillators: 



X„ 



-V'(X n ) + C{X n ^+X n 



2X„ 



(19) 



where the oscillator potential V possesses one minimum, 
V'(0) = and V"(0) = 1. The spectrum of small ampli- 
tude plane waves is given by w 2 = 1 + AC sin 2 (q/2) and 
discrete breather solutions are time-periodic spatially lo- 
calized solutions of Eq. (19) with frequency fib ^ uo q /m 
for any nonzero integer m. Here we will again consider 
a weak interaction between the oscillators (C « 1), 
and assume that the breather is essentially a single-site 
excitation X (t) = X (t + 2n/fl b ) ^ and X n ^ {t) « 0. 



The equations for the linearized phase space flow around 
the breather solution are then given by [6] 



e'n — -£n + C(e n -i + e n+ i — 2e n ) 
-6 n>0 (v"(X (t)) - l) e • 



(20) 



For C = we expand V"(X (t)) = J2k v k e 
the Floquet representation e (t) = J2k e ofe e 
arrive at the set of equations, 



fc ° b * and use 

i(ui+kQ b )t j- 



-(lu + kfl b ) 2 e k 



E 



(21) 



for the site n — 0, at which the breather is excited. This 
complete set of linear equations describes the coupling of 
the open channel (the corresponding amplitude of oscil- 
lations eoo) and many closed channels (the amplitudes of 
oscillations e fe with k ^ 0). 

To obtain the condition for total reflection we apply a 
procedure, which is similar to our analysis of the DNLS 
system. We "turn off" the coupling between the open 
channel and closed channels, i. e. the amplitude eoo is 
set to zero in Eq. (21). The remaining homogeneous 
set of equations determines the local states wl of closed 
channels. The values of u = ujl arc those for which the 
corresponding determinant vanishes. Moreover, we are 
interested in the situation when w| = 1, which results 
in a resonance of a local mode of the system of closed 
channels with the open channel oj 2 w 1. It is evident 
that such a situation is not necessarily realized for an 
arbitrary value of fib. To find the proper value of fib we 
put uj — 1 and scan fib for a given potential V (note 
that the parameters Vk depend on fib). As a result, we 
expect to find a discrete set of fib values, for which the 
determinant of the reduced set of equations Eq. (21) 
zeros, and a Fano resonance occurs. 

Thus at variance to the DNLS Klein- Gordon 

chain in the anticontinous limit will provide us with a 
total reflection of waves by breathers only for a selected 
discrete set of breather frequencies. Increasing the cou- 
pling C transforms each of these frequency values into 
frequency stripes on the real axis, which will continue 
to increase with C. This follows from the fact that the 
band width of uo q increases linearly with the coupling C, 
whereas the shift of the eigenvalues is proportional to C 2 . 

Note here, that the Green function technique elabo- 
rated for the wave scattering by DBs in Ref. [6], leads 
to the same procedure described above in the limit of 
small Vk, allowing to obtain the dependence of ujl on the 
amplitude and frequency of the DB. 

In order to demonstrate the validity of our constructive 
approach, we choose 



V(x) 



1 



1 



1 



(22) 



and carry out the above procedure. We obtain fib « 1.38. 
The prediction then is that for small C a breather with fib 
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close to this value yields a total reflection. The numerical 
result for the transmission is shown in Fig 2 for C — 0.001 
[14]. We indeed observe a total reflection around q « 1, 
as predicted. 



energy (E F = E n ). 
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FIG. 2. Transmission coefficient 

versus q for a Klein-Gordon chain breather with Qb = 1.38 



and C = 0.001. 

Our analysis can be directly used for the detection of 
Fano resonances in the propagation of linear waves in 
weakly dissipative Josephson junction arrays in the pres- 
ence of discrete breathers [2], [6]. Moreover our results 
allow to formulate optimal conditions for electronic trans- 
port spectroscopy of molecules, quantum dots or optical 
cavities, coupled to a system of leads. Indeed a one 
dimensional electronic transport through the localized 
states in the presence of externally applied microwave 
radiation of frequency f2, will display a Fano resonance 
(and correspondingly zero conductivity) as the condition 
Ep = fi ± E n is satisfied. Here, Ep is the Fermi energy, 
and E n are the localized levels. To establish this type 
of a spectroscopy leads should be quasi-one-dimensional 
such that only one transversal mode persists. Then we 
reduce the lead system to one open channel. The other 
important condition is to satisfy that the wavelength of 
excitations (e.g. at the Fermi energy for electronic trans- 
port) is large compared to the spatial extension of the 
coupling to the localized states. Additional optional spa- 
tial modulations, e.g. in the leads, may be used to gener- 
ate artificial gaps in the electronic spectrum and thus to 
tune the wavelength. The power of microwave radiation 
has to be rather small to exclude the renormalization of 
E n . The coupling between propagating states and lo- 
calized states may be also provided by, e.g. spin-spin or 
spin-orbit interactions. In this case zero conductance will 
be the consequence of a bare dot state passing the Fermi 
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